Quantum-disordered ground state for hard-core bosons on the frustrated square lattice 
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We investigate the phase diagram of hard-core bosons on a square lattice with competing inter- 
actions. The hard-core bosons can also be represented by spin-1/2 operators and the model can 
therefore be mapped onto an anisotropic Ji- J2-Heisenberg model. We find the Neel state and a 
collinear antiferromagnetic state as classical ordered phases to be suppressed by the introduction 
of ferromagnetic exchange terms in the x-y plane which result in a ferromagnetic phase for large 
interactions. For an intermediate regime, the emergence of new quantum states like valence bond 
crystals or super-solids is predicted for similar models. We do not observe any signal for long-range 
order in terms of conventional order or dimer correlations in our model and find an exponential 
decay in the spin correlations. Hence, all evidence is pointing towards a quantum-disordered ground 
state for a small region in the phase diagram. 

PACS numbers: quantum spin frustration, 75.10.Jm; Magnetic phase transitions, 75.30.Kz; computer mod- 
eling and simulation 75.40.Mg 



I. INTRODUCTION 

The investigation of frustrated spin models has become 
a rather active field over the past years due to the rising 
interest in now quantum phases like valence bond solids 
or spin liquids. ^'^ Another motivation to study these sys- 
tems is the question of the microscopic origin of high- Tc 
superconductivity which is still controversally discussed 
and often connected with frustrating spin interactions. 
One of the most interesting and challenging problems 
in this field is the J1-J2 -spin- 1/2 isotropic Heisenberg 
model. ^'-^"^ There is still no final answer to the question 
of the intermediate phase in the ground state phase dia- 
gram of this model and several techniques have been used 
to track it down. Non-variational Quantum Monte-Carlo 
(QMC) simulations have a severe sign problem for the 
frustrated model and are, hence, very limited for such a 
system. In this paper we present our work on an anal- 
ogous bosonic model which maps onto the Heisenberg 
model^*^"^^ for a certain set of parameters and may give 
some crucial hints for the completely frustrated model. 

Starting from the classical model without quantum 
fluctuations, which was analyzed in the early eighties by 
Landau and Binder^^^^^ and was further investigated in 
recent years, ^^^^^ we examine the quantum model for fi- 
nite temperatures and extrapolate our QMC results to 
T = to draw a ground state phase diagram. 

For the equivalent anisotropic Heisenberg model we 
find two classical magnetically ordered phases (Neel and 
collinear state) as ground states for small quantum fluc- 
tuations and a direct transition between these two anti- 
ferromagnetic configurations. For large fluctuations the 
system becomes ferromagnetic in the x-y-plane. Close 
to the highly frustrated point which is accompanied by 
a large ground-state degeneracy in the classical limit we 
flnd a region with no finite order parameter and inter- 
pret this state as quantum-disordered. In the bosonic 
language the antiferromagnetic states are described by 
boson-density waves with wave vectors q ~ (tt, tt) and 



q = (0, tt) or (tt, 0) respectively. The ferromagnetic in- 
plane order is interpreted as Bose condensation of the 
magnons and, hence, corresponds to a superfluid order in 
the bosonic model. ^^"^^ From now on we will use these 
terms exchangeably. 

The paper is structured as followed: in the subsequent 
part we introduce the model and give an overview of the 
underlying physics and critical points of the system. In 
the third section we discuss some insights on the model 
via perturbation theory. The fourth section is divided 
into three subsections and dedicated to the QMC simu- 
lations. In the first part we briefiy introduce the QMC 
method we used and explain the difficulties with the ther- 
malization process for the frustrated model. Thereafter 
we present and discuss our observations of the magnetic 
and quantum-correlated observables that yield the phase 
diagrams at finite and zero temperature. In section five 
wc show results from an exact diagonalization for one set 
of parameters. In the concluding part of the paper we 
will discuss our results and give an outlook for further 
calculations. 



II. MODEL 

The Hamiltonian of the model is given by summation 
over all interactions of nearest neighbors (NN) and next- 
nearest neighbors (NNN) and the according exchange in- 
tegrals for these bonds: 
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The are bosonic creation (annihilation) operators 

with (6-^'')^ = (hard-core bosons) and rn = 6|6j is the 
occupation number for the site i (limited to or 1). The 
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model resides on a N = L x L square lattice with peri- 
odic boundary conditions (for a sketch see Fig. 7 below). 
The Vi > are chosen to be repulsive and the hopping 
integrals ti < are negative. Thus, mapping the bosonic 
operators onto spin-1/2 operators, the model is equiva- 
lent to the anisotropic Ji-J2-spin-l/2 Heisenberg model 

NN NN 

+ J2 {s+sr + h.c.) + j| s^s^ (2) 

NNN NNN 

using ti = J^'^ /2 and V = Jf. This model is frus- 
trated in the component (antiferromagnetic interac- 
tions, Jf > 0) and not frustrated in the and 5^ 
components (ferromagnetic interactions, Jf'^ < 0). For 
2ti = J^'^ = Jf = Vi the isotropic Heisenberg model is 
recovered. Since we are only interested in the case of a 
half-filled bosonic model {N/2 bosons for the whole lat- 
tice) we do not take into account a chemical potential in 
(1) and, hence, the magnetic field in the corresponding 
antiferromagnetic Heisenberg model (2) is zero and we 
are working in the subspace of (5*^) = 0. Earlier works 
calculating the bosonic model with ti < - as we do 
here - have focused on the effect of a varying chemical 
potential. ^""^'^ Furthermore, Refs. 10-12 did not consider 
next-nearest neighbor hopping. 

In the classical limit, i.e., for ti 0, the Hamiltonian 
represents the antiferromagnetic Ji- J2-Ising model which 
is highly frustrated in the region J2 « j^/2.^^'^^'^°'^^'^^ 
The phase diagram contains two magnetically ordered 
ground states - Neel order for J2 < Ji/2 and coUinear 
order for J2 > Ji/2 - and the paramagnetic phase for 
high temperatures. At J2 = Ji/2 the ground state is de- 
generate (of order 2^+^) which leads to a suppression of 
the critical temperature and freezing problems in simple 
Monte-Carlo simulations. For small quantum fiuctua- 
tions I til we expect the classical ordered phases to sur- 
vive for low temperatures and to build the ground state. 
At the frustrated point V2 = Vi/2 we calculated second 
order perturbation terms by means of degenerate per- 
turbation theory (DPT) to classify the behavior in the 
vicinity of the critical point (see part HI). For large \ti\ 
we expect a superfluid phase which corresponds to long- 
range ferromagnetic order in the x-y plane in the spin 
model. 

For the most interesting region (intermediate \ti\ and 
V2 ~ Vi/2) the emergence of new quantum states 
was predicted by Balents et al?'' for similar models. 
The smallest building block for these quantum states 
is a dimerized configuration of two spins which in our 
anisotropic model needs a short explanation. In the ba- 
sis of total spin 5**°* of these two spins there exist two 
entangled states with {S^) = given by: 

i5o) = = 0, = 0) = (in) - i ;t)) , 
\s{) = |5*°* = 1, = 0) = 4 (in) + 1 ;t)) . (3) 




Figure 1. (Color online) Ground state phase diagram cal- 
culated with QMC simulations. As abscissa we plotted the 
degree of frustration V2/V1 and as ordinate the magnitude of 
the quantum fluctuations —ti/Vi (with ti/V\ = t2/V2)- The 
dashed blue hue indicates the direct transition from Neel to 
colhnear order derived by means of DPT (up to second or- 
der, see section III). The gray area represents roughly the 
region where we do not find any finite signal for various order 
parameters. 



The energies for these two eigcnstates depend explicitly 
on the amplitude t of the quantum fluctuations: 

{t{S+Sj + h.c.) + VStS^^)\S^) = (-t - y/4)|5o) , 
{t{StS- + h.c.) + VSfS^)\S{) ^{t- y/4)|5i) . (4) 

Hence, for our anisotropic model with ti < 0, the low- 
est energy is given by the ^i-dimer state. The covering 
of the whole lattice with dimers allows for different or- 
dered configurations as, e.g., a columnar or staggered 
valence bond solid or for non-static arrangements as a 
resonating valence bond solid. Apart from these dimer- 
ized states a superposition of the classical ordered states 
- order in and in-plane order at the same time - is 
feasible. In the bosonic language these states are called 
supersolids.^^'^'^'^^'^^"'^^ 

As the main result which will be derived below we al- 
ready show the ground state phase diagram in Fig. 1 
where the suppression of the two magnetically ordered 
phases can be seen. Furthermore, it can be seen that the 
onset of the ferromagnetic phase for large fluctuations 
does not coincide with the opening of the two antiferro- 
magnetic phases. Hence, we find an intermediate phase 
without conventional order. Furthermore, we present the 
finite-temperature phase diagrams for different magni- 
tudes of quantum fluctuations in Fig. 2. It can be seen 
that the classical ordering process is suppressed to lower 
temperatures due to the kinetic energy introduced by the 
exchange integrals ti. 
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only second order corrections to the classical energies and 
using t2 ~ ti/2 and V2 ~ Vi/2 - yields the relation: 

1^.1 [iV2 3 

which is shown in the final ground state phase diagram 
(Fig. 1) as dashed line for small Since equation 

(5) does not depend on the sign of U it holds also for 
antiferromagnetic x-y interactions and can be compared 
to the result of a series expansion by Oitmaa et al.^^ 
where for small fluctuations ti > the direct transition 
between the classical antiferromagnetic states survives as 
well. 



Figure 2. (Color online) Phase diagrams at finite tempera- 
tures for different values of ti. For increasing \ti\ the transi- 
tion temperatures are shifted to lower values and conventional 
order is suppressed. As a guide to the eye the full (blue) line 
shows the known phase boundary for the classical appropi- 
ately scaled Ising model. ^° For ti = —0.1 Vi (broadly dashed 
(green) line) only transition temperatures close to the critical 
point were calculated. 



III. DEGENERATE PERTURBATION THEORY 

To estimate the influence of small quantum fluctua- 
tions on the classical ground states nearby the critical 
point we calculated the second-order perturbation in the 
degenerate ground state manifold at the critical point. '^^ 
For V2 = Vi/2 every state with 2 bosons per plaquette 
(square of 4 lattice sites) has the same classical energy. 
Thus, for the whole lattice which is constituted of over- 
lapping plaquettes this local degeneracy yields a global 
one of the order of the lattice length L, as explained 
in reference 20. The DPT distinguishes between diago- 
nal perturbations which leave the system in exactly the 
same state and off-diagonal perturbations which transfer 
the system into another state of the degenerate mani- 
fold. In the case of the frustrated square lattice different 
ground states are connected via flips of antiferromagnetic 
spins in a whole line or row of the lattice. Thus, the or- 
der of off-diagonal perturbations scales with the length 
of the lattice L. Off-diagonal perturbations are therefore 
negligible. However, diagonal perturbations are already 
relevant in second order and non-constant for the two 
classical starting points - Neel and coUinear configura- 
tions. In the Neel state only nearest-neighbor hopping 
ti is possible on 2 bonds of the lattice and yields an 
energy gain of AEi = —2 L^r^. In the collinear state 
nearest-neighbor hopping on bonds and next-nearest- 
neighbor hopping on 2 bonds is possible which gives an 

^2 ^2 

energy gain AE2 = —2 i^jpj — L'^^. Thus, small fluc- 
tuations enforce the classical ground states in the vicinity 
of the critical point. Calculating the transition line be- 
tween the Neel and collinear state - taking into account 



IV. QUANTUM MONTE-CARLO 
A. Algorithm 

For the calculation of the complete phase diagram at 
finite temperatures and in the ground state we use quan- 
tum Monte-Carlo techniques. For negative hopping inte- 
grals ti the Stochastic Series Expansion (SSE)'^'^"'^^ has 
no sign problem and the directed loop algorithm^^ yields 
an adapted update scheme for a large set of parameters 
in this model. We used an implementation of the ALPS- 
project'^^"^^ as basis for our simulations. 

However, the frustration in the model produces a crit- 
ical slowing down in the Monte-Carlo (MC) simulation 
and the large degeneracy in the vicinity of the critical 
point V2 = ^1/2 causes severe thermalization problems in 
the standard implementation. To overcome these prob- 
lems we added an exchange MC update^*^^^^ in tempera- 
ture space to the ALPS directed-loop apphcation. Start- 
ing from the expansion of the partition function in the 
SSE scheme one can derive an appropriate Metropolis up- 
date probability due to detailed balance for this update. 
To guarantee a good thermalization within the exchange 
MC step it is important to adjust the temperature steps 
and number of sweeps between the exchanges of config- 
urations (swaps) of neighboring simulations. 

In addition to the exchange MC, wc used an anneal- 
ing procedure for each copy independently during the 
thermalization process. This kind of algorithm helps to 
prethermalize the simulations at lower temperatures to 
ensure a better swap rate for the exchange MC algorithm. 



B. Magnetic order 

To determine the regions in the phase diagram where 
the system is classically ordered, we calculated the struc- 
ture factors 
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Figure 3. (Color online) Temperature dependence of the 
fourth-order cumulants for the Neel order parameter for dif- 
ferent lattice sizes at V2 = 0.2 Vi and ti — —0.15 Vi. In the 
inset the estimate of the error for the transition temperature 
is given. 
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Figure 4. (Color online) Comparison of the superfluid density 
(spin stiffness) ps (right) and cumulants for magnetic order 
parameters (left) - Neel order (NO) and coUinear order (CO) 
- for —ti/Vi = 0.5 and V2/V1 — 0.5. The clear signature in 
the superfluid order parameter indicates a finite-temperature 
phase transition. 



for Neel order {q = (tt, tt)) and coUinear order {q = (tt, 0) 
and (0, tt)) for various parameters of the Hamiltonian (1) 
to check for antiferromagnetic order in the S"^ direction. 
To detect the exact transition temperature into the mag- 
netic phases we took the fourth-order (Binder) cumulants 
of the order parameter:^'^'''^ 

which intersect for different lattice sizes L at the critical 
temperature. As an example the scenario at V2 = 0.2 Vi 
and ti = —0.15 Vi is shown in Fig. 3. Using this method 
we constructed the finite-temperature phase-diagrams for 
different ratios of quantum fluctuations —ti/Vi for the 
two antiferromagnetic phases and the high-temperature 
paramagnetic phase shown in Fig. 2. The transition 
temperatures are suppressed to lower values for increas- 
ing hopping terms. For —ti/Vi > 0.1 and values of 
V2 ~ Vi/2, i.e., in the vicinity of the critical point no 
magnetic ordering in the 5*^ direction can be detected. 

We extrapolated the finite-temperature results to T = 
to draw the phase boundaries of the magnetic phases 
in the ground state phase diagram in Fig. 1. The shape 
of the phase diagram is very similar to that calculated 
by Oitmaa et al.^"^ for the anisotropic Ji- J2-Heisenberg 
model with antiferromagentic interactions in S^'"^'^ by 
means of series expansion in ti around the Ising limit. 

Furthermore, we performed measurements for the su- 
perfluid density (or in spin language the spin stiffness)^"* 
which is measured via the variance for the net-direction 
of the off-diagonal operators in the QMC simulations.^^ 
This quantity indicates a correlated movement of the 
hard-core bosons or a ferromagnetic order in the x-y 
plane respectively. For large S^, S'^ interactions \ti\ we 
find this state to be the ground state. For —ti/Vi = 0.5 




Figure 5. (Color online) Ground state values of the or- 
der parameters for Neel order (NO), coUinear order (CO), 
in-plane ferromagnetic order (FO) and the energy density 
at ti = —0.25 Vi and varying frustration V2/V1. (Sym- 
bols are larger than error bars and values do not change for 
20 > L > 16). 



and V2 = 0.5 Vi the temperature dependence of ps is 
shown in Fig. 4 for two different lattice sizes L = 8, 14. 

In Fig. 5 we show the evolution of the order pa- 
rameters (Neel, coUinear and in-plane order) for differ- 
ent frustration strengths and fixed quantum fiuctuations 
—ti/Vi = 0.25 at a sufficiently low temperature. The 
measurements are converged to their ground state values 
and do not depend on the lattice size any more. The tran- 
sitions from Neel order to in-plane ferromagnetic order 
and back to antiferromagnetic order (coUinear configura- 
tion) are clearly visible and we conclude from the sharp 
features in the order parameters that the transitions are 
of first order. However, calculating the observables for 
smaller values of \ti\/Vi gives rise to interesting behav- 
ior especially for the ferromagnetic order parameter. For 
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Figure 6. (Color online) Top: Ground-state values of the 
order parameters for Neel order (NO), coUinear order (CO), 
in-plane ferromagnetic order (FO) and the energy density at 
ti = —0.1 Vi and varying frustration V2/V1 (symbols are larger 
than error-bars). Bottom: Evolution of the ferromagnetic 
order parameter ps at the critical point V2 = Vi/2 for in- 
creasing ratios \ti\/Vi. For small fluctuations ps ~ even 
after the direct transition between Neel and collinear order is 
suppressed; only for |ti|/Vi > 0.175 a finite value of the order 
parameter is measurable. (Symbols are larger than error bars 
(if not given) and values do not change for 20 > L > 16) 



the highly frustrated region around V2 ~ Vi/2 and small 
\U\/V^ (with values 0.08 < \U\/Vi < 0.175) the signature 
of Ps depends on the lattice size and goes to zero for low 
temperatures and large lattices. Thus, we find a small 
region without any of the conventional order parameters 
giving a clear sign for an ordering process (Fig. 6 top). 
In addition we show the trend of the ferromagnetic order 
parameter at V2 = Vi/2 in Fig. 6 (bottom) where the cal- 
culations in the intermediate region needed larger lattices 
and very low temperatures (20 > L > 12, T « 0.01 Vi) 
to converge to their ground state value. In smaller lat- 
tices a strong tendency to superfluid order was noticed 
due to the periodic boundary conditions. For increasing 
lattice size the order parameter began to oscillate be- 
fore vanishing completely for small temperatures. Since 
this calculation is very time-consuming we performed it 
only exemplary for the single value of V2 = Vi /2 and do 
not give the exact phase transition from the disordered 
phase into the ferromagnetic phase in Fig. 1 for all values 
of ¥2- However, the smooth behavior of the order param- 
eter and the energy in Fig. 6 indicates that the transition 
into the ferromagnetic phase is probably continuous. 




Figure 7. Local measurements of magnetization and kinetic 
bond energy in a grayscale for a 12 x 12 lattice at V2 = 0.51 Vi 
and ti = -0.1 Vi in the ground state (T = 0.01 Vi). No 
ordered features are distinguishable. 



C. Quantum correlations 

The lack of the expected conventional order in a fi- 
nite region of the ground state phase diagram motivated 
further simulations and calculations of new order param- 
eters. We performed measurements for the whole struc- 
ture factor and thereby ruled out any magnetic order in 
the S'^-direction. We also checked for unconventional or- 
der, i. e., quantum-ordered phases like columnar or stag- 
gered dimer phases. For this purpose we implemented 
the measurements for bond-bond correlations which cor- 
respond to four-spin correlation functions. In the SSE 
one can use an improved estimator for correlations be- 
tween operators that are part of the hamiltonian itself. '^^ 

The measurements of local quantities as local magne- 
tization and local kinetic bond energies are shown in 
Fig. 7 for parameters which lie in the critical region 
{~ti/Vi — 0.1 and V2 = 0.51 Vi) in a grayscale. There 
is no distinguishable magnetic order and the energies of 
the bonds are equally distributed. Diagonal bonds have 
smaller energies due to the fact that \t2\ < \ti\. 

The correlation measurements of spins and bonds are 
shown in Fig. 8. As a reference we chose the top left spin 
with its right adjacent bond. All other bonds and sites 
represent the strength of the correlation from this special 
site (bond) to the reference site (bond). For the bond- 
bond correlations kinetic and potential energy terms are 
taken into account and they are normalized to the auto- 
correlation of the reference bond which yields the largest 
value (black in the figure). There is no long-range or- 
der to identify neither in spin-spin correlation which de- 
cay exponentially fast nor in the bond-bond correlation 
which decay all equally to a finite value that is given by 
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Figure 8. We show correlations between spins and dimers liv- 
ing on nearest-neighbor bonds on a periodic 12 x 12 lattice, 
hence up to 6 neighboring spins and bonds. We chose the top 
left spin and adjacent right bond as reference points for the 
measurement. All other sites and bonds of the plot represent 
the correlation of spins and dimers in the illustrated distance 
to the reference bond coded in a grayscale. For V2 = 0.51 Vi 
and ti = —0.1 Vi and in the ground state (T = 0.01 Vi), i.e., 
in the region without conventional order, spin correlations de- 
cay exponentially to zero and dimer correlations decay equally 
to a finite value, hence, we find no long-range order neither 
in dimer nor in spin correlations. 



the local bond energy. Here the weakest correlations are 
measured for orthogonal bonds of the same square (top 
left) where the strongest correlation is given by the par- 
allel bond on the same square. We conclude from these 
measurements that there does not exist any long-range 
order in the system for a finite region of parameters. 



V. EXACT DIAGONALIZATION 

A topologically ordered state could be another pos- 
sibility for the phase without any signatures for order. 
To check for this kind of non-local ordering, the calcula- 
tion of the spectrum is necessary since a degeneracy of 
the groundstate is expected for a topologically ordered 
state. ^'"^^ The spectrum is not accessible via QMC simu- 
lations and therefore we performed an exact diagonaliza- 
tion for a 6 x 6 lattice with periodic boundary conditions 
- as on a torus - at tj = —0.1 Vi and V2 = 0.5 Vi. 
For the computation of the spectrum we used an exist- 
ing implementation of an exact diagonalization by Jorg 
Schulenburg.^'' 

We calculated the lowest eigenvalues in the different 
/c-sectors and show the energy difference to the lowest 
eigenvalue Eo{0, 0) = —10.2854 Vi in Fig. 9. For symme- 
try reasons, it is sufficient to concentrate on the region 
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Figure 9. Exact diagonalization data for a 6 x 6 lattice with 
periodic boundary conditions at V2 ~ 0.5 Vi and ti = —0.1 Vi. 
Shown are the energy differences AE = Ek — Eo for differ- 
ent fc-values (only ky < fca,) in a grayscale (AE'min(7r, tt) — 
0.1609 Vi). The spectrum shows no tendency for degener- 
acy and therefore a topologically ordered state is unlikely. In 
the inset the same values are shown along the path (0, 0) — )■ 
(7r,0) (7r,7r) (0,0). 



< kx < TT, < ky < kx of the Brillouin zone. The 
main panel of Fig. 9 shows a grayscale plot of the energy 
differences in this region, the inset a different representa- 
tion of essentially the same data. We obtain a minimal 
energy gap AE^^in = 0.1609 Vi for k = (tTjTt) above a 
unique ground state on the 6x6 lattice. This is compa- 
rable to the dispersion along the fc^-dircction where we 
find a maximum gap AE{2tt/3,0) = 0.3482 Vi. Accord- 
ingly, we interpret the large value of Ai^niin as evidence 
against a ground-state degeneracy. 

The spectrum shown in the inset in Fig. 9 is also qual- 
itatively different from the spectrum of the Heisenberg 
model (cf., e.g., Rcfs. 4 and 48). Indeed, in the latter 
case one would expect a behavior close to fc^, = 0, fcj, = 
which is similar to the vicinity of the ordering wave vec- 
tor (which is kr^ = TT, ky = TT for the Neel state). The 
excitation spectrum shown in Fig. 9 is therefore consis- 
tent with the absence of classical order as demonstrated 
by the QMC simulations. 

While it would be desirable to perform a finite-size 
analysis, we have selected the 6x6 lattice as the only 
accessible lattice which has the symmetries of the infi- 
nite system and is compatible with the expected ordered 
states (the 4x4 'square lattice' is not necessarily rep- 
resentative for two dimensions since it can also be inter- 
preted as a four-dimensional torus). Still, the gap Ai^min 
is comparable to the dispersion of the excited states even 
on the 6x6 lattice which suggests that the gap will also 
stay finite in the thermodynamic limit. According to this 
result no groundstate degeneracy and for this reason no 
topological order is present. 



7 



VI. DISCUSSION 

We simulated a two-dimensional lattice model for 
hard-core bosons with competing interactions which 
maps onto an anisotropic frustrated quantum spin model. 
We used improved QMC methods to calculate classical 
magnetic order parameters for a large parameter space 
of the model. We found two antiferromagnetic phases 
as expected from the classical limit (the frustrated Ising 
model) and a ferromagnetic configuration which corre- 
sponds to a superfluid phase in the bosonic language. 
Most importantly, we detected a finite region without any 
magnetic order. Careful calculations of higher correla- 
tions as four-spin correlations to check for dimer ordered 
phases which were proposed by Balents et al. for similar 
models^^ did not give any signal of long-ranged order. 
Further exact diagonalization calculations showed that 
the system is not topologically ordered, either. In con- 
clusion, there seems to exist a finite region in parameter 
space without any long-range order and with an exponen- 
tial decay in the spin-spin correlations which has not been 



observed in earlier works on similar models. -'^'^^^'^ Since we 
did not find any dimcr-ordered phases we assume that 
their appearence strongly depends on the kinetic energy 
encoded in the model. An introduction of different ex- 
change terms as, e. g., a ring exchange could enhance the 
probability of finding valence bond crystals. ^''''^^ 
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